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The field and temperature dependence of the de Haas- van Alphen oscillations spectrum is studied 
for an ideal two-dimensional compensated metal. It is shown that the chemical potential oscillations, 
involved in the frequency combinations observed in the case of uncompensated orbits, are strongly 
damped and can even be suppressed when the effective masses of the electron- and hole-type or- 
bits are the same. When magnetic breakdown between bands occurs, this damping is even more 
pronounced and the Lifshits-Kosevich formalism accounts for the data in a wide field range. 
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I. INTRODUCTION 

In large enough magnetic field, the Fermi surface (FS) 

I of multiband quasi- two-dimensional metals, is liable 
to give rise to networks of orbits coupled by magnetic 
breakdown (MB). The most studied type of network is 
the linear chain of coupled orbits introduced by Pippard^ 
and illustrated by several quasi-two-dimensional (q-2D) 

! organic conductors such as k-(BEDT-TTF)2Cu(NCS)2. 
As discussed, in Refj^, magnetic oscillations spectra 

\ of such networks contain many frequencies that are 
linear combinations of two basic frequencies. In ad- 

\ dition to those linked to MB-induced orbits, other 
frequencies are observed that are not accounted for by 
the semi-classical theory of Falicov and Stachowiak^. 
They can be attributed to quantum interference (as 
far as magnetoresistance oscillations are concerned), 

' MB-induced modulation of the density of states^i^i^ and 
oscillation of the chemical potentialii^iiSiiiS, even though 

\ the actual respective contribution of these three phenom- 
ena to the oscillatory behavior remains to be established. 

Another type of network is provided by q-2D metals of 
which the FS is composed of compensated electron- and 
hole- types tubes. This is the case of the family of or- 
■ ganic metals (BEDT-TTF)8Hg4Cli2(C6H5X)2 (X = CI, 
\ Br) whose FS, which originates from two pairs of cross- 
ing q-lD sheets, is composed of one electron and one 
' hole tube with the same areaii. As it is the case of the 
above-mentioned linear chains of coupled orbits, magne- 
toresistance oscillations spectra in this type of network 
reveal frequencies that are linear combinations of three 
basic frequencies, linked to the compensated orbits and 
to the two FS pieces located in-betwee n^^i^'^i-'^^ . How- 
ever, in striking contrast to the data relevant to linear 
chains of orbits, de Haas- van Alphen (dHvA) oscillations 
spectra recorded in the case of the compound with X 
= Br only exhibit oscillations, the field and tempera- 
ture dependence of which can be consistently interpreted 
on the basis of the semiclassical model of Falicov and 



Stachowiak^ii^. Analogous networks are observed in or- 
ganic metals, with two carriers per unit cell. In this case, 
the FS originates from the overlapping in two directions 
of hole tubes with an area equal to that of the First Bril- 
louin zone (FBZ) and from the resulting gap openingsi^. 
As reported in the case of (BEDO)4Ni(CN)4-4CH3C]Nfi^, 
such a FS yields a network consisting in two hole- and 
one electron-type tubes (see Fig. [DJa)). Closely related 
network is obtained in the case where the large hole or- 
bit come close to the FBZ boundary, as it is the case of 
(BEDT-TTF)4NH4[Fe(C2 04)3]^. In this latter case, a 
large MB gap is observed at this point and the resulting 
network only consists in one electron- and one hole-type 
orbit, as displayed in Fig. [TJb). Linear chain of successive 
electron-hole tubes might also be observed in Bechgaard 
salt (TMTSF)2N03i^iia. 

The aim of this paper is to explore the field and tem- 
perature dependence of the dHvA oscillations spectra of 
an ideal 2D metal whose FS is composed of one electron- 
and one hole- type compensated orbit. It is demonstrated 
that the field-induced oscillations of the chemical poten- 
tial are strongly damped for such a FS and can even 
be suppressed in the case where electron- and hole-type 
orbits have the same effective masses. The chemical po- 
tential oscillations can be even more damped when the 
two orbits are coupled by MB. In this case, the Lifshits- 
Kosevich (LK) formalism accounts for the data up to 
high magnetic field and low temperature. 



II. MODEL 

We first consider a 2D metal whose electronic structure 
consists of two parabolic bands of hole- and electron- type, 
respectively. The bottom of the electron band is set at 
zero energy while the top of the hole band is at A > 0. 
The total number of electrons in the system is such that 
the hole band is completely filled. Since the lower part of 
the electron band is lower in energy than the top of the 
hole band, some quasiparticles move to the electron band 
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FIG. 1; (color online) Calculated Fermi surface (FS) of q- 
2D organic metals (a) (BEDO)4Ni(CN)4-4CH3CNi& and (b) 
(BEDT-TTF)4NH4[Fe(C2 04)3]^ leading to networks of com- 
pensated electron- and hole-orbits (solid red and blue lines, 
respectively). Ellipses in dotted blue lines correspond to the 
hole orbits, with an area equal to that of the FBZ, from which 
the FS is built (see text). 



in order to lower the total energy. The effective masses 
linked to the two bands, m* = 1/Ce and m\ = 1/Ch > 0, 
can be different. It is useful to define the physical units 
of the problem. The reduced field b = eBAo/h is in 
units of the characteristic field B ~ h/eAo, the effective 
masses are in units of the electron mass rrie, the energies 
are in units of _B = 2'!Th^ /uieAo and the temperature t in 
units of T = E/ks- Given a unit cell area =197.6 A^, 
which stands for (BEDT-TTF)8Hg4Cli2(C6H5Cl)2^, we 
obtain 13 = 2093 T and f = 2812 K. Therefore, real- 
istic experimental conditions yield small values of h and 
t compared to B and T, respectively. The semi-classical 
quantization of the energy levels in the presence of a mag- 
netic field leads to the Landau equations: 



Ee{n) 



Ceb{n+-), 



Eh{n) 



Chb{n+-) (1) 



each Landau level (LL) having a degeneracy h per sam- 
ple area. The zero field Fermi energy is simply given 
by Ep = rn^l (m* + m^JA. At finite temperature, the 
total free energy is given by the difference between the 
contribution of the electron and of the hole band, with 
the condition that the total number of quasiparticles 
Npj-i is fixed. In addition, the compensation condition 
imposes that the number of quasiparticles in the elec- 
tron {Ne) and in the hole {Nh) band arc the same. 
From the thermodynamical relations, we can define a 
free energy for the system based on the difference be- 
tween the free energy for the electrons in the electron 
band and the free energy for the holes in the hole band 
/S.F = fie — Oft, + [Ne — Nh)iJ,, where fle{h) is the Grand 
Potential for the electrons (holes) and /i is the chemical 
potential [^{t = 0, 6) = Epib)]: 



n,{t,b) = -tbJ2^og{l + cxp[f3{^i-E,{nm 

ri>0 

nhit,b) = tbJ2^og{l + exp[P{Ehin)-fi)]) (2) 



Since = Nh, we conclude that F — ^l^ — ^h hi com- 
pensated metals, fi is evaluated from the self-consistent 
equation dAF/dn = 0. At zero temperature, the above 
expressions reduce simply to the ground state (GS) en- 
ergy AEq. The Fermi energy Ep{b) is given by the con- 
dition rieb = rihb, where rig and are the (integer) num- 
bers of LL filled, and 



Tle-l 



AEo ^ -bJ2 iEpib) ~Ee{n))-bY, (Ehin) - Epib)) 



E„ 



Eh + b{nh - n,)EF{b) = E,^ - Eh 



(3) 



The special cases where the Fermi energy goes trough 
one Landau level, or where this Landau level is partially 
filled, correspond to singular points in the energy spec- 
trum as a function of the inverse field that do not modify 
the thermodynamical quantities. The exact expression 
for A_Eo is simply 



AE^^^{Ce + Ch)b''nl 



bAuf, 



(4) 



with Uf, = [Ep/CJ} + l/2]i, the notation standing 
for the integer part of the argument. The GS energy 
oscillates around the limit of zero field, where Ai?o(<7 = 
1) = -A2/2(Ce -f- C/i). For example, taking A = 1, 
Ce = 1 and Ch = 2/5, we obtain AEq ~ -0.357. We 
deduce the oscillating part of the magnetization Mosc — 
—dAEo/db from the latter expression, using succesively 
the Fourier transforms of the periodic functions [x]i — x 
and [[x\i — xY: 



= -(Ce + Cft)Fo V sin(2^fci^o/&) (5) 



k>l 



where Fo = A/(Ce + Cft) 



i,:A/(TO: 



I) is the 



fundamental frequency corresponding to the FS area of 
the electron and hole band. At zero temperature and for 
fixed number of electrons, the magnetization oscillates 
with characteristic frequency Fq, and the amplitudes Ak 
of the A:*'' harmonics are given by the LK formula with 
1/k depcndcnco2i: 



Ai 



-Ch 



(6) 



The sum + Ch means that the 2 orbits circulating 
around the hole and electron bands contribute individu- 
ally to the magnetization. 
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III. SELF-CONSISTENT EQUATION FOR THE 
CHEMICAL POTENTIAL 



The oscillatory parts of the Grand Potentials Eqs. ^ 
can be extracted using Poisson's formula for any function 
F{x) (see Ref,^^ for details): 



1 r°° 

FiO) + Y,F{n)^ / Fix)dx 

^ n>l 

25RV / F{x)exp{2iTrnx)dx 
„ ^ , Jo 



(7) 



n>l ' 



where F{n) is equal to log(l + exp[/3(/i — Ee{n))]) or 
log(l + exp[f3{Eh{n) — fi)]) for electrons and holes, re- 
spectively. The last series in Eq. ([7]) gives the oscil- 
latory part of the Grand Potential in terms of Fourier 
components. The index n is expressed as a function of 
the energy n = n{E) by relations |T]), then expanded 
around the chemical potential ^ at low temperature 
where the energy derivatives of the distribution functions 
1/(1 + cM(i{Ee{n) - or 1/(1 + exp[/3(Ai - -B/.(n))]) 
are strongly peaked. After some algebra, we obtain for 
each band 



1 

2 
1 



(8) 



X °° f— 1)" u 



12 ^ TT-n 



(9) 



+ E ^^("^D cos(27rn— 
12 TT^n^ Uhb 

n—l 



where 



-R("™e(/i)) 



(10) 



is the temperature reduction factor for effective 
masses nm*^^j^y A Dingle term RD{nm*f^f^yt*f^f^j) = 
exp(— 27r^nm*^^^t*^^^/&) can also be added in the case 
where the relaxation times t* and tjj for electron- and 
hole-band have to be taken into account. In the follow- 
ing we will assume that these two relaxation times are 
negligible for convenient purpose. It is always possible 
to add those terms in final expressions. The chemical 
potential satisfies therefore the self-consistent relation: 
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FIG. 2: (color online) Field dependence of the chemical po- 
tential normalized to its value in zero-field at t = 0.001. Solid 
and dashed lines correspond to one electron and one hole com- 
pensated orbits (with m^ = 1 and = 1.1 or 5/2) and two 
electron orbits (with mS — 1 and mi — 1.1 or 5/2), respec- 
tively (see text). 



At a first order approximation, which will be discussed 
in the following sections, we can replace /x in the sine 
functions of the previous expression by Ep when the os- 
cillations of the chemical potential are small compared to 
Ep^ in the small field and high temperature regime {t/b 
large). This gives: 



Ef 



ml 



ml 



(12) 



Y^^—J-[R[nml) 

n=l 



■KTl 



i?(nm*)] sin(27rn 



It can be remarked first that, in the case where the ef- 
fective masses and relaxation times linked to the electron- 
and hole-type orbits are the same, there is an exact solu- 
tion for Eq. pT|) with /i = Ep . In this case the oscillatory 
part of Eq. (fTTj) or Eq. (fT2|) vanishes and the chemical 
potential remains constant in magnetic field and temper- 
ature. This is due to the fact that the energy levels are 
symmetric around Ep. More generally, the amplitude of 
the chemical potential oscillations can be compared to 
the case of two electron bandsi^. In Fig. [51 the field- 
dependent chemical potential Ec^. pT|) is calculated for 
compensated orbits with m* = 1 and mj^ = 1.1 or 5/2 
and compared to the case of two electronic orbits with 



effective masses mn = 1 and 



1.1 or 5/2. It can 



be observed that the chemical potential oscillations are 
strongly damped for compensated orbits, even in the case 
where m^ and m* have strongly different values^^. 



E 



(-1)" A 
—[R{nm*i^) sin(27r?T, 



-h 



n=l 



i?(nm*) sin(27m—)] 



Chh ' 
(11) 



IV. DE HAAS- VAN ALPHEN OSCILLATIONS 

The oscillatory magnetization can be obtained putting 
the solution of the chemical potential given by Eq. (fTTj) 
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F/F„ 

FIG. 3: (color online) Fourier spectrum of the magnetization 
at t = 0.001 for m* = 1 and m^ = 1.1 or 5/2. The field 
range is from b = 0.015 to b = 0.033. Fq is the fundamental 
frequency (see text). 



back into the expression for the free energy F using Eqs. 
^ and As discussed above, the chemical potential 
is field-independent for nie = m.h and the LK formula3i 
holds in this case. However, for m* ^ m'^, the oscillatory 
magnetization differs from the LK theory. Examples 
of Fourier analysis of the magnetization are given in 
Fig. [3] It should be noted that, contrary to the case of 
electron-hole systems away from exact compensation^, 
frequency combinations are not observed directly in 
the oscillatory spectra. However, as discussed later on, 
deviations from the LK behavior are observed for the 
second harmonics. 




It is useful for experimentalists, to check to what extent 
it is possible to determine a temperature and field range 
in which the LK formalism provides a satisfactory ap- 
proximation of the oscillatory behavior in the case where 
m* 7^ m^ . The b/t dependence of the Fourier compo- 
nents of the magnetization with frequencies Fq and 2Fo 
are given in Fig. |4l The LK formula accounts for the field 
and temperature dependence of the first harmonics. Fur- 
thermore, it can be noticed that in the case where mj^ is 
strongly different from m*, the contribution of the orbit 
with the smallest effective mass dominates in a large field 
range (see the dashed line in Fig. 0^). It has also been 
checked that, in the opposite case where m* is close to 
m^ , the data can be accounted for by the contribution of 
only one orbit with a mean effective mass, namely, m*„g^„ 
~ (m* -I- m*f^)/2. In contrast, the LK formula cannot ac- 
count for the second harmonics in the case where m* and 
m^ significantly differ. The observed behavior of the sec- 
ond harmonics in this latter case is reminiscent to that 
observed for two electronic orbitsi^. 



V. MAGNETIC BREAKDOWN 

In this section, we consider the case where the mag- 
netic field is large enough for the quasiparticles can tun- 



FIG. 4: t/b dependence of the Fourier components of the 
magnetization with frequencies Fo and 2Fo. The effective 
mass of the electronic-type orbit is ml = 1. The effective 
mass of the hole- type orbit is m^ = 5/2 and 1.1 in (a) and (b), 
respectively. Solid lines are obtained from the LK formula. 
The dashed line in (a) corresponds to the contribution of the 
electronic orbit, only. 



nel through MB between the two bands. 

The topology of the new FS is depicted in Fig. [S] 
the quasiparticles can tunnel between the junctions a' 
and P (/?' and a) with probability amplitude ip or be 
reflected between the junctions a' and a (/?' and /3 re- 
spectively) with probability amplitude q. This surface is 
a simplification of Fig. [TJ) in the sense that it does not 
include the ID network feature, p and q are related by 
p^ + q^ = 1——, and the field dependence of p is given by 
the Chambers formula p = exp(— where bo 
is a characteristic MB field, bo comes in the field range 
covered by experiments in the case where the two bands 
are located closely enough in the FBZ. Depending on 
the band type (electron or hole), the quasiparticle wave- 
function acquires a phase either (jg = Se/b or cr^ = Sh/b 
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FIG. 5: Fermi surface of a two-band model with magnetic 
breakdown in the FBZ. The quasiparticles orbit clockwise for 
the Hole surface and counterclockwise for the Electron band. 
They can tunnel between the points (which represent also the 
amplitudes of the wave-function) a' and /3 (or and a) with 
a probability amplitude ip, or be relected between a' and a 
(or /3' and /9) with a probability amplitude q. At the turning 
points (cross symbols), the wavefunction acquires a factor i 
(or —i) depending wether the quasiparticle orbits an electron 
or hole surface. 

when orbiting around each part of the Fermi surface (be- 
tween a and a' counterclockwise for the electron band or 
f} and (3' clockwise for the hole band). The scmiclassical 
actions Se = 2Trm*E and Sh = 27rmJJ (A — E) represent 
the areas delimited by each Fermi surface in the FBZ. 
Assuming that the wavefunction of the quasiparticle is a 
at the point a, we obtain that after one orbit around the 
electron band, a' = — aexp(i(Te) (the quasiparticle goes 
through 2 turning points where it acquires a factor i or 
— i each time, depending on the electron or hole typeS^). 
Similarly, it follows = — /3exp(— icr/i) for the hole band. 
The four amplitudes a, a', (3 and /3' are related by the 
transfer matris^i^ 

This set of equations have non-zero solution when the 
energy E satisfies the implicit equation: 

(1 + qexpiae){l + qexp{-iah)) 

+p^expi{ae-ah) ^0 (14) 

In the case where the bands arc disconnected {p — 
0), Eq. can be factorized and reduces simply to 

1 -I- exp(i(Te) = or 1 -I- exp(— i(T;i) = 0, which corre- 
sponds to the two independent sets of discrete energy 
levels given by Eq. ([1]). In general Eq. can be 

solved numerically for p ^ 0: the field dependence of the 
LLs energy is plotted in Figs. [H] and [7] for q=Q.6 and 
for 2 different sets of effective masses. In the case where 
p ^ 0, gaps open at the former LL intersections, and the 
structure of the two bands is modified. Alike in the case 
where p = 0, a quasi-hole and a quasi-electron band can 
be defined. Indeed, increasing the magnetic field, levels 
that go upwards above the Fermi energy Ep (b) (see dia- 
mond symbol lines) and downwards below Ep(h) repre- 
sent quasi-electron and quasi-hole bands, respectively. At 
zero temperature, the compensation condition N^, ~ 
is replaced by the condition that the lower quasi-hole 
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FIG. 6: (color online) Energy levels for a reflection amplitude 
q = 0.6 (or, equivalently, tunnelling amplitude p — y^l — q'^), 
A = 1, m^ = 1 and mj = 5/2. Levels that belong to a quasi- 
electron band (going upwards in the limit of large field, black 
lines) above the Fermi energy Epip) (diamond symbols, green 
line) can be distinguished from levels that belong to a quasi- 
hole band (respectively going downwards, red lines). The 
Fermi energy at zero field is shown as a dashed line. 

band is filled completely, the upper electron-band being 
empty. The two bands are therefore always separated by 
the chemical potential energy, as shown in Figs. [5] and 
[71 As the magnetic field decreases, the field-dependent 
Landau levels near the Fermi energy become flatter since 
the levels of each band do not intersect the Fermi level. 
In the case where q — Q [p — 1)^ the quasiparticles tunnel 
directly through the 2 bands at each revolution, and the 
exact solution of Eq. ([14]) is given by the set of LLs: 

i?e,.(<7 = 0,n) = Ep-b-^^^{n + \) 

with n > 0. Within the semi-classical approach, the 
cross-section area of the FS corresponding to this in- 
verted parabolic quasi-hole band, which is completely 
filled, is zero. Since it is proportional to the frequency of 
the oscillations, there are no oscillations at g = 0. This 
simple feature is in line with the prediction of the semi- 
classical model for two compensated orbits, keeping in 
mind that the area of these orbits has opposite signs^. 

We assume that Cu = CePa / qa where po a-nd go sxc co- 
prime intcgers^iiS. Indeed, it is useful, also for computing 
time reasons, to approximate any real value of Ch/Ce by 
a close rational ratio since in this case the spectrum 
is periodic in energy, with periodicity Tp = Cebpo 

In order to estimate the total free energy in the gen- 
eral case, it is necessary to solve numerically Eq. [M] We 
assume that Ch = CePo/qo where po and go arc coprime 
integers^iiS. Indeed, it is useful, also for computing time 
reasons, to approximate any real value of Ch/Ce by a 
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FIG. 7: (color online) Energy levels for a reflection amplitude 
q = 0.6, A = 1, m* = m'^ = 1. Levels that belong to a quasi- 
electron band above the Fermi energy located at Ep{h)=Q.^ 
(diamond symbols, green line) can be distinguished from lev- 
els that belong to a quasi- hole band (red lines) . The 2 bands 
are fully symmetric with respect with the chemical potential 
energy 



close rational ratio since in this case the spectrum (|14p is 
periodic in energy, with periodicity Te ~ Cebpo Within 
each interval Te, there are exactly po + go discrete so- 
lutions (there is indeed conservation of the number of 
LLs when q decreases from unity). By successive energy 
translations, it is then easy to reproduce all the spec- 
trum. The GS energy AEq given by Eq. ^ is however 
no more well defined, due to the mixing between electron 
and hole levels. Indeed, the GS energy should corre- 
spond to the sum of all the quasi-hole energies below the 
Fermi energy Epib). These energies being not bounded 
by below, the GS energy is formally infinite. We propose 
to account numerically for this problem by introducing a 
cutoff function in the density of states with the necessary 
condition that the magnetization should not depend on 
variations of the cutoff parameters. 

Given a set of LL energies Eehici,n-), n = 0, • • ■ ,po + 
go — 1, solutions of (|14[) in any energy interval of width 
Te, we can introduce a cutoff function ifdE) such as 
(pc{E) = 1 for £■ larger than a characteristic energy Ec 
and equal to exp[— c(i? — Ec)^^] for E < Ec, where S is 
any positive integer greater than 1 (we take S ~ 4: in the 
simulations which gives a very smooth cutoff function). 
This function has the required property of conserving 
the important physical properties near the Fermi surface 
and making in particular the GS energy finite. The LL 
density of states Pc{E) takes the following form 



Neh = / dEpciE) 



rEpib) Po+go-1 oo 

= b I dE V ME) 



dE 

5{E- E,h{q,n)-kTE) 



n— A.;— — oo 



(16) 



where N^h is, as mentioned before, the total number 
of quasiparticles in the Canonical Ensemble. Once the 
parameter c is obtained, we can define for example the 
GS energy Ai?o 

/Epib) PQ+qa-l oo 
dE 
-oo n 



n— k— — oo 

E5{E - Ech{q,n) ~ kTE) 



(17) 



or the free energy AF 



/oo Po-l-90-1 oo 
-oo n I 



3-1-90-1 oo 

Ai^ = -tb 

'-00 n=0 k=-c^ 

log[l + expf3{p - E)]S{E - Eeh{q,n) - kTE) 
Neh4^ (18) 



The potential p{t, b) is calculated from Eq. ([T8|) by ex- 
tremizing the free energy dAF/dp = 0, and compelling 
the magnetization Mosc — —dAF/db to be independent 
of the parameter Ec for Ec far away from the chemi- 
cal potential or at energies large compared to the Lan- 
dau gap. We have checked, for different values of Ec 
in units of A, for example, Ec = —1,-2, —4, and for 
a large range of fields, that the resulting magnetization 
does not change. We choose Neh, which is arbitrary, as 
a multiple of the characteristic zero field energy density 
m* -I- = (Ce + Ch)/CcCh times the zero field Fermi 
energy Ep. For Ec = —2, we take in the following nu- 
merical simulations Neh = SEpiCe + Ch) I CeCt- Then, 
for each value of the field b, the parameter c defined from 
Eq. (|16p is unique. Also, in the case q ~ 1, the numerical 
solution for the magnetization is fully consistent with the 
results of the first section in absence of magnetic break- 
down. In Fig. [51 the oscillations of the chemical potential 
fi{t, b) obviously decrease with q until the total tunneling 
occurs where only orbits with frequency zero are allowed. 



A. Analytical amplitudes for small field 



Po+qa-i oo 

Pc{E) = b Y ME)5{E-Ech{q,n)-kTE) 

n=0 fe= — oo 

Given Ec, the parameter c is found to be solution of 
the equation of conservation at zero temperature: 



In this section, we extract the analytical expression 
for the first harmonics amplitude (corresponding to the 
frequency Fq) in the small b/t regime, and compare it 
to the numerical results of the previous section. We first 
assume that the oscillations can be described by an effec- 
tive free energy which is constructed by adding reflection 
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FIG. 8: (color online) Oscillations of the chemical potential 
But t = 10~^ for different values of g; A = 1, = 1 and 
= 5/2. 
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FIG. 9: (color online) Comparison between the numerical so- 
lution for /i (Eq. I18|l with approximation (a) from Eq. (|20|l . 
and (b) from Eq. (|20|l with jj, in the sine functions replaced 
by Ef- ml = 1 and = 5/2, and q = 0.6. The 2 sets 
t — 0.004 and t = 0.008 are shifted along the vertical axis for 
clarity. 



amplitudes to the temperature reduction factors in 
the expression of the Grand Potentials Eqs. [5] and [51 for 
orbits that circulate n times around Fermi surfaces Se 
or Sii- In this approximation, it is indeed clear that the 
quasi-particles orbit the same surface if the field is small 
enough, and do not tunnel through the junction. The ef- 
fective free energy can then be derived directly from Eqs. 
dSl) and (El) as 



AFeff = Qe- 



(A - 



2Ce 2Ch 



(19) 



b^Ce ^ (-g)" 



— 7T-^R{nm*,) cos(27rri,— — . 



2 ^ — ' TT^n 

n—l 



and the chemical potential is derived as in Eq. [11] 



fi ^ Ef + V [R{nml) sin(2^n ^ ^ 

^ 71 — 1 



R{nml) sin(27rn-4^) 



Chb ' 
(20) 



Numerically, we can measure the deviations between 
this approximation and the numerical result [i ~ 
dAF/db, where AF is defined by Eq. dH]), for m* = 1 
and ml = 5/2 (see Fig. [9]), and find that Eq. ^ is 
valid for b/t smaller than approximatively 12 (approxi- 
mation (a)), and 8 in the case where is replaced by Ep 
in the right hand side of Eq. (|20l) (approximation (b)), 
as in Eq. HH) for q=l. 

The complete derivation of the first harmonic ampli- 
tude is done in Appendix A. If we keep the dominant 
terms in Eq. (|A7p . which are those with the dominant 



reduction factor Ri'Til^h))^ ^^'^ discard the other expo- 
nentially small ones R{nm*^^^f^-^) for n > 2 or products 
"^e(/i))^' obtain 



where a, 



TT \ m* Ue ml ah J 

eih) - 2ml^^^q{R{m*)~R{ml))/{ml+ml) as 
defined in Appendix A, and Ji is the Bessel function of 
order 1. In the limit where ae(h-) £^rc small, the functions 
2Ji(Q;e(^))/ae(/i) are very close to unity and we recover 
the LK formula: 



LK 



Fo(R{ml) , i?(m;j 



(22) 



For example, for q = 0.6, to* = 1 and to^ = 5/2, a, 
vanishes when b/t is either large or small, a, being al- 
ways less than 0.2 elsewhere. At this extremum value, 
the function 2J(x)/x is approximatively equal to 0.995. 
Therefore we do not expect, as in the case g = 1, a sig- 
nificant deviation from the LK formula at small b/t, as it 
can be seen in Fig. [10] in this range of fields for different 
values of q strongly different from unity. 



B. LK formula for large field 

In the large b/t limit, the previous approximation is no 
more valid (the reduction factors are close to unity). We 
show hereafter that there are an infinite number of orbits 
contributing to the frequency Fq. These orbits, which are 
described in details in Appendix B, have large effective 
masses me(n) = {n+l)ml+nml or mh(n) = nm^ + {n + 
1)to^, where n is a positive integer. Their contribution to 
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the total amplitude are therefore negligible at low field 
(or high temperature) since their reduction factors are 
exponentially small. An example of orbit of mass TOe(l) 
contributing to frequency Fq is shown in figure 1111 

We can however compute the exact amplitude Ai 
within the LK theory, and show that there is no dis- 
cernible difference with the numerical solution of Eq. 
in the Canonical Ensemble. 

At low fields, the first harmonics amplitude Ai con- 
tains the contribution of the 2 single orbits of mass m* 
and mj^, as in Eq. ([22]) . At higher fields, we can replace 
the previous quantity by an exact expression 

AL1< = EliT^^Th) (23) 

TT 

where the amplitudes re(/i) includes all the orbits that 
contribute to the first harmonics 

Te{h) = -|— i?(m:(^)) 

^1 TOe(/i)(n) 

The combinatorial coefficients ie(/t) ('^) defining the to- 
tal amplitude for nonequivalent orbits of mass ?7ig(/j) (n) 
are computed in Appendix B. By definition, we set 
te(h) (0) = q- For comparison, we have solved numeri- 
cally the magnetization from Eq. for different values 
of q and extracted the first harmonics Ai from Fourier 
analysis. In Fig. [TO] the amplitude of the first harmon- 
ics Ai is plotted versus b/t for q = 0.6,0.8,0.98 and 1 
(symbols). Since q depends on h through the Chambers 
formula, these plots should merely be regarded as the 
temperature dependence of Ai at a given magnetic field 
value. The data are compared to the predictions of the 
LK theory (sohd lines) given by Eqs. ^ and (|B5|) . 

for m* = 1 and = 5/2. The deviations with the LK 
formula are negligible in all the h/t range explored, what- 
ever the q value. In particular, the contributions of the 
higher mass orbits which become important in the large 
field range are still well described by the LK formula. 



0.3 



~ 0.2 

CL 



. 1 i i i 1 i i 1 1 1 i i i i 1 1 . 

^ - • — — • •— 




• — • — 


\q = 0.98 


/ rf^*** — 


= 0.8 




y — q 


= 0.6 




■ . . . . 1 . . 


. 1 


. . . . 1 . 



50 100 150 

b/t 

FIG. 10: (color online) First harmonics amplitude A\ for dif- 
ferent values of g; A = 1, m*^ = 1 and = 5/2. The filled 
circle symbols represent the numerical analysis from the free 
energy expression psp . and the lines represent the LK for- 
mula (|23l I24|l . with elements te(h)('i) of Eq. (|B10p computed 
up to n = 10 from Eq. (|B5|) . 




FIG. 11: Example of orbit contributing to first harmonics Fq. 
The eff'ective mass here is me(l) = 2m* -f mj^. This orbit can 
be represented by the operator QePP or PPQe (see Appendix 
B for details). The point Oe on surface Se is crossed twice by 
the trajectory. 



VI. SUMMARY AND CONCLUSION 

The field-dependent chemical potential oscillations in 
FS composed of two compensated electron and hole orbits 
is strongly damped when compared to the case of a FS 
with only electronic orbits^iSiiSiiiS. It is even suppressed 
in the case where the effective mass of electron (m* ) and 
hole (mjjj) band are the same (assuming the relaxation 
times t*i^ii-f are either identical or negligibly small). In 
addition, the LK formula accounts for the field and tem- 
perature dependence of the first harmonic's dHvA oscil- 
lations amplitude in all cases. As for the amplitude of the 
second harmonics, it is accounted for by the LK formula. 



provided m* and m^ are not strongly different. Besides, 
as MB develops, the chemical potential oscillations are 
further damped and the previous conclusion on the LK 
validity still holds for the amplitude of the first harmonic. 
In this well as in the general case of compensated 

electron-hole orbits, the contributions of an infinite num- 
ber of orbits to the main frequency become relevant at 
high fields, due to the existence of closed trajectories in 
the FBZ which have zero frequency and make the ex- 
pression of the LK amplitude rather complex. Within 
the FS topology that we have considered here, we have 
demonstrated that it can nevertheless be computed ex- 
actly. We expect these additional " zero frequency" orbits 
to be reminiscent in other systems as well. 
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Finally, it should be noticed that the Landau energy 
spectrum of ID and 2D periodic networks of electron-hole 
compensated orbits (see Figs. [T]^b) and (a), respectively), 
is likely a set of Landau bands rather than discrete Lan- 
dau levels as observed in Figs. [6] and [T] This feature 
could induce a MB-induced modulation of the density 
of states, as observed in non-compensated networks^^i^ 
and hence frequency combinations in dHvA oscillatory 
spectra. In that respect, the reported method for solving 
Eq. [13] can be used for the study of such orbits networks. 



For extracting the first amplitude, we select the 
integers such as to + 71 = ±1 in (jA4p and to — 
n = ±1 in (|A5[) . In the product G{x)G'{x) of ex- 
pression (jA2p . the terms can be rearranged noticing 
that cos(27rTO^0''c) sm{2'!TnFox) = [sin(27r(n + m)Fox) + 
sin(27r(n — m)Fox)]/2 with m,n > 1. For the first har- 
monic, only terms n — m = ±1 contribute. A further 
approximation is to keep the first term of the series (jAip . 
i.e. G{x) w -q{R{ml) - R{ml)) sm{2TT Fqx) / 1: , in ([M]) . 
so that the coefficients i?e(/i) {"^^ can be computed ex- 
actly, from the explicit relation 



APPENDIX A 



From Eq. (|20p . wc can define the periodic function 



exp [2m'u;g(/i)g(i?(m*) — i?(TO^)) sin(27rFoa;)] = 

00 

Jm(nae(/i))cxp(2i7rTOFoa;) (A6) 



m— — 00 



G{x) = - — '^^[R(nml) — R(nml)] sin(27rni^o2;l(Al) where J™ is the Bessel function of order m and a, 



where x ~ 1/b. The chemical potential Eq. ([20]) be- 
comes /i = Ep + bG{x)/{m* + TO^). Replacing /i in 
pop by the previous expression, we compute the oscil- 
lating part of the magnetization AIosc = —dAFef f /db ~ 
x^dAFeff/dx: 



G{x)G'{x) 



CeCh 
~ [Ce + Ch) 

1 V i-iT 

n=l 

ChR(nmt)-^^exp(2iTrn ^ ^ ^ x) 
ox Ch 



(A2) 



Cg i?(nm* ) — K exp (2z7rn — x) 
ox Og 



We also define the Fourier coeffieicnts BQf}^\{n^m) as 



+00 



2i7TmF()X 



(A3) 



where We(^h) = C'/i(e)/(Ce + Ch)- Then, we perform the 
derivatives in expression (jA2| by noticing for example 
that 



— sRexp(2i7rn-^x) = 3? 

ox Ce 



——00 

X 2i7r(m + n)Fo exp(2i7r(m + n)Fox)] (A4) 



+00 



and 



— 5Rexp(2i7rn^^^— x) = K 
ox Ch 



Bhiji.m) 

.— — 00 

X 2z7r(m — n)Fo exp(2z7r(m — n)F^x)\ (A5) 



+00 



e(/i) - 

2we(/j)(7(i?(TO*)— i?(TO^)). In this case it is easy to identify 
Be{h){n,rn) = Jm(nae(ft))- After some algebra, we find 
that the amplitude Ai for the first harmonic Fq is given 
by the expression, in the large t/b limit 



■kAi 



V2i::^((-i)"c.i?(nTO:):^^ 



n=l 



^ r>/ *-,Jninah)\ , CeCh 

ChR{nml) — ^ > + — — X 



nah 



Ce+Ch 

Y,-[R{nml)^R{nm*^]x (A7) 



{[i?((n-l)TO,:)-i?((n-l)TO:)]g2"-i 
-[i?((n + 1)to;;) - R{{n + l)m:)]q^-+'] 



APPENDIX B 

In this appendix, we compute the amplitude of the first 
hamonics in the LK theory for any given reflection am- 
plitude q. In the lowest approximation, the contribution 
to the amplitude is given by one single orbit around the 
electron or hole band (see Fig. O. Each contributes with 
a mass equal to to* = l/Cg and to^ ~ 1/Ch respectively 
and a reflection amplitude q at the junction point. At 
finite temperature, we add a reduction factor R[m*) or 
R{rri^) so that we obtain (with a Fq factor overall) Eq. 

However, there are other contributions to the Fq har- 
monics, coming from more complex orbits. For example 
in Fig. II1[ from a starting point Og near the jimction on 
the Fermi surface Se, we can reach the junction point 
and go through to the other band St, then complete 
an entire orbit and come back through the same junc- 
tion to the band Se, and finally perform 2 orbits around 
until we reach again. The frequency will be pro- 
portional to 25*6 — Sh- For a compensated system, this 
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frequency is just Fq since Se = Sh- The total mass is 
however proportional to the derivative of 2S'e — with 
respect to the energy, which is 2m* + m'^ (the derivative 
of Sh = 2TTm'^{A — E) is indeed negative). Therefore we 
should take into account all possible trajectories. This 
can be done exactly by summing up all the contributions 
for a given mass me(n) = (n + 1)to* + nmjj in > 0). As 
before we start from a point Oe (Oh) on the surface Se 
{Sh) and turn until we reach the junction. Here wc note 
Qe (Qh) the reflection operator which keep the quasipar- 
ticlc on the surface Se with amplitude q (Sh respectively), 
and P the operator which take the quasiparticle through 
the junction (with amplitude ip). Their electron-hole 
representations are given by the matrices: 

In the simplest case, the orbit is described only by the 
operator or graph Qe, whereas in the second example, 
the orbit is described by the graph PPQe- The orbit 
QePP is equivalent and has the same mass 2m* + mj^. 
Indeed, the orbit goes trough the same point Oe twice, 
and therefore there are 2 different ways of writing the 
graph above, depending on which branch of the orbit we 
place the point Oe- Equivalent graphs are described by 
a cyclic permutation of their operator elements. We also 
observe that multiple orbits (PP)" have zero frequency 
and mass n(m* +m^J. We will note in the following te{n.) 
{th{n)) the sums of the amplitudes of all the trajectories 
with the same mass ■me{n) (mh{n) = nm* + (n + i)m'^ 
respectively), which are constructed starting from an ar- 
bitrary point Oe (or Oh respectively). Te (Th) will be the 
total amplitude for the frequency Fq {—Fq) starting from 
point Oe (Oh) on the surface Se (Sh), see Eqs. (|23I24|) 

Th is computed in the same way as Te by replacing 
the label e in all the quantities by h. For equal masses 
m* = mj^, we have clearly Te = Th- For example, the 
orbits and amplitudes contributing to the mass me(l) = 
2m* -|- mj^ are 

PPQe ^ QePP ^ (V)?, 

This gives te{l) = {—p^)q- The 2 operators above cor- 
respond to the same graph by a cyclic permutation of 
their elements PP and Qe and the sign ^ is understood 
as an equivalence between identical graphs. In the next 
case, for the mass me (2) = 3m* -I- 2ml^ we have the pos- 
sibilities 

PQhPQl ~ QePQhPQe - QlPQhP ^ ( 

PPPpQe - PPQePP - QePPPP ^ {-p'^fq 

we obtain te{2) = {—p^q^ +p'^q). In the general case, 
we can construct all the possible orbits that contribute to 
the mass me{n) by representing them as multiple prod- 
ucts of elementary operators 



Q''e'PQ'f^'PQZ^PQ'^^'-P--Q"e''PQh''P (B2) 

where > and m^ > are positive integers. The 
constraints imposed by the mass on theses integers and 
k are 

k k 

Ui = n — k + 1, TTii = n — k, 1 < k < n (B3) 

1=1 1=1 

There are also k cyclic permutations of the graph (jB2[) 
above, corresponding to moving elements Qe'PQ^h'^ 
the right or to the left. For example, an equivalent graph 
would be 

QT-PQrP--QTPQh''PQ7PQrP (B4) 

The total amplitude te {n) for a given n can be written 

as 

te{n) = f^c„,fe(V)'=<z'^"-'^+' (B5) 

fc=i 

where the coefficients Cn,k enumerate the number of 
all nonequivalent graphs having the same mass, up to 
cyclic permutations. It corresponds to all possible sets of 
TLi and rrii representing graphs (jB2[) with the constraint 
(|B3|) . divided by the number of cyclic permutations k. 
The total amplitude te(n) can be written more precisely 
as 

n 

te{n) ^ J2 J2 qSti«.+i:Li™.(_p2)fc X 

k—l rii ,mi >0 

'k^l2'l=i n,,n-fc+l^E?=i m,,n-k (B6) 

The Kronecker functions Sk,o can be represented by 

integrals 6kfi — Jq cxp{2i'Kkx)dx — <f dz/{2iTrz) x z'^, 
and the last expression becomes 




Summing up the series, we obtain 



" (-p2)fc f dz z'^-"-^ r .^lk-n-1 

~ k^i / 2i^ (1 - qzY f 2A^{1- qz'Y 

Applying the residue theorem, it is easy to show that 
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and 



dz 



2i7r (1 — qz)'' 



\k-ii 



l\k ^ \k-l)1 



2in (1 - qz') 



(B8) 



(B9) 



where (jj) = n\/kl{n — fc)! are the binomial coefRcients. 
We obtain finally 



TABLE I: Values of the coefficients c„_fc which represent 
the number of non-equivalent orbits for a given mass (n + 
l)m*(^) + ^i^h(e) with 2k breakdowns, 1 < A; < ti. 



n \ 


1 


2 


3 


4 


5 


6 


7 


Mass 


1 


1 
















2 


1 


1 












3ni*(h) + 2m^(e) 


3 


1 


3 


1 










4ni*(h) + 3m^(g) 


4 


1 


6 


6 


1 








S'Ti^fh) + 4m*(^) 


5 


1 


10 


20 


10 


1 






6»7i*(h) + 5m,;(^) 


6 


1 


15 


50 


50 


15 


1 




7ni*(fe) + 6m;;(^) 


7 


1 


21 


105 


175 


105 


21 


1 


SnT-:^,,, + 7m*f^, 



Cn,k 



■n \ fn-l\ 
.k-l) \k-l) 



(BIO) 



The values of the coefficients c„_fc up to n = 7 are given 
in Table 1. 



We also observe by symmetry that th{n) is equal to 
te{n) since this geometric coefficient does not depend ex- 
plicitely on masses. 
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